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Abstract. We present a phenomenological jet model which 
is able to reproduce well the observed variations of the 
submillimetre-to-radio emission of the bright quasar 3C 273 
during the last 20 years. It is a generalization of the original 
shock model of Marscher & Gear (1985), which is now able 
to describe an accelerating or decelerating shock wave, in a 
curved, non-conical and non-adiabatic jet. The model defines 
the properties of a synchrotron outburst which is expected to 
be emitted by the jet material in a small region just behind 
the shock front. By a proper parameterization of the average 
outburst's evolution and of the peculiarities of individual out- 
bursts, we are able to decompose simultaneously thirteen long- 
term light-curves of 3C 273 in a series of seventeen distinct out- 
bursts. It is the first time that a model is so closely confronted 
to the long-term multi-wavelength variability properties of a 
quasar. 

The ability of the model to reproduce the very different 
shapes of the submillimetre-to-radio light curves of 3C 273 
gives strong support to the shock model of Marscher & Gear 
(1985). Indirectly, it also reinforces the idea that the outbursts 
seen in the light-curves are physically linked to the distinct fea- 
tures observed to move along the jet with apparently superlu- 
minal velocities. The more than 5000 submillimetre-to-radio 
observations in the different light-curves are able to constrain 
the physical properties of the jet. The results suggest, for in- 
stance, that the magnetic field behind the shock front is rather 
turbulent. There is also some evidence that the jet radius does 
not increase linearly with distance down the jet or, alternatively, 
that the synchrotron emitting material decelerates with distance 
and/or bends away from the line-of-sight. 

Key words: radiation mechanisms: non-thermal - galaxies: ac- 
tive - galaxies: jets - quasars: individual: 3C 273 - radio con- 
tinuum: galaxies 



1. Introduction 

The theory of synchrotron emission in relativistic jets was de- 
veloped around 1980 mainly by Blandford & Konigl (1979), 



Marscher ( 198C ) and Konigl (1981 ). In 1985, both Marscher & 



Gear (1985 hereafter MG85) and Hughes et al. (1985) pro 



posed a model for the emission of a shock wave propagat- 
ing down a simple relativistic jet. Both groups assumed that 
the jet is confined to a cone of constant opening angle and 
that t he jet fl ow is adiabatic. The computer code of Hughes 
et al. ( 1989a ) was able to describe the low-frequency flux den- 
sity an d polarization variability of BL Lacertae (Hughes et al. 



1989b). The model of MG85 has however the advantage to in- 



clude the effects of synchrotron and inverse-Compton energy 
losses of the electrons, which cannot be neglected at higher 
frequencies. It provides a simple explanation for the distinct 
components in the jet observed using very long baseline inter- 
ferometry (VLBI). Multi-wavelength total flux measurements 
were found to be difficult to use to constrain the shock model 
of MG85, because the emission of all distinct features in the jet 
often overlap to for m a nearly flat total spectrum, as illustrated 
by Marscher ( |1988[ ) for the quasar NRAO 140. 

It is only since 1995 that the very well sampled total flux 
millimetre and radio light-curves of a few sources allowed 
to study the spectral evolution of individual synchrotron out- 
bursts. By subtracting a quiescent spectrum assumed to be con- 
stant, Litchfield et al. (|1995|) and Stevens et al. (\199% 11996, 



1998) could follow the early evolution of single synchrotron 
outbursts in 3C 279, PKS 0420-014, 3C 345 and 3C 273, re- 
spectively. These studies gave additional support to the MG85 
model, but failed to constrain it strongly, mainly because they 
could only follow the evolution of an outburst until the onset of 
the next one. 



To overcome this problem, Tiirler et al. (1999t, hereafter 
Paper I) proposed a new approach to derive the observed prop- 
erties of synchrotron outbursts, which consists in a complete 
decomposition of long-term multi -wavelength light-curves into 
a series of self-similar flaring events. Two different approaches 
are presented in Paper I to describe the evolution of these events 
with both time and frequency. The so called "light-curve ap- 
proach" is model-independent and describes empirically the 
shape of the light-curves of individual outbursts at different 
frequencies. A second approach based on three-stage shock 



models, like those of MG85 or Valtaoja et al. (1992), describes 
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directly the evolution of the flaring synchrotron spectrum. Fi- 
nally, the "hybrid approach" of Tiirler (2000) models strictly 
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the shape of the synchrotron spectrum, but leaves the evolution 
of the spectral turnover as free as possible. The results of this 
third approach, based on self-similar outbursts as in Paper I, 
were found to be in very good agreement with the spectral evo- 
lution expected by the shock model of MG85. 

After this last test, the next step, which is presented here, is 
to adapt the shock model of MG85 in order to describe physi- 
cally both the average evolution of the outbursts and their indi- 
vidual specificity. This generalized shock model, described in 
Sect. 1^, takes into account the effects of an accelerating or de- 
celerating synchrotron source in a curved, non-conical and non- 
adiabatic jet. This model is confronted to the observations by 
fitting seventeen distinct outbursts simultaneously to thirteen 
submillimetre-to-radio light-curves of 3C 273. The proper pa- 
rameterization to achieve this fit is described in Sect. ^ and the 
results of this light-curve decomposition are given in Sect. ||. 
In Sect. ^ we discuss the implications of our results, both on 
the global properties of the inner jet and on the peculiarities of 
individual outbursts. The main results of this study are summa- 
rized in Sect. ^ 

Throughout this paper, we use the convention Si, oc for 
the spectral index a and we use "Ig", rather than "log", for the 
decimal logarithm "logj^Q", because of a lack of space in tables 
and long equations. 



these fluxes as a free parameter of our fit.[] To smooth out the 
dips in the GBI light-curves (cf. Tiirler et al. 1999ai, we aver- 



age all GBI observations into bins of 10 days. Finally, we do 
not consider the GBI measurements before 1980, to avoid a flux 
increase in the light-curve due to an outburst that started before 
1979. We end up with a total of 5234 observational points to 
constrain the shock model. 



3. Physical shock model 

In the original shock model of MG85, the shock front was as- 
sumed to propagate with constant speed in a straight conical 
and adiabatic jet. Partial generalizations of this model were al- 
ready derived by Marscher (1990) and Marscher et al. (1992) 
in the case of a bending jet, and by Stevens et al. (1996) in 
the two cases of a straight non-adiabatic jet and a curved adia- 
batic jet. Here we further generalize the shock model of MG85 
to account for the effects of an accelerating or decelerating 
shock front in a curved, non-conical and non-adiabatic jet. In 
Sect. [3.1[ we describe the typical three-stage evolution of all 
outbursts, whereas in Sect. 3.2 we show how the physical con- 
ditions at the onset of the shock can influence the evolution of 
individual outbursts. 



2. Observational material 

The light-curves fitted here are part of the m ulti-wa velength 
database of 3C 273 presented by Tiirler et al. ( 1999a ). The 12 
light-curves from 5 GHz to 0.35 mm are as described in Paper I, 
except that we now consider the observations up to 1999, in- 
cluding the most recent measurements from the Metsahovi Ra- 
dio Observatory in Finland and from the University of Michi- 
gan Radio Astronomy Observatory (UMRAO). 

We extend to lower frequencies the analysis of Paper I 
by adding a new light-curve at 2.7 GHz. This light-curve is 
constituted of observations from the Green Bank Interferom- 
eter (GBI) and from the 100 m telescope at Effelsberg in Ger- 
many (Reich et al. 1998). The fluxes obtained with the GBI 
at 2.7 GHz and 2.25 GHz include both 3C 273B (flie inner jet) 
and 3C 273A (the hot spot at the far end of the jet), but these 
two components combine partially out of phase (Ghigo F.D., 
private communication). As a consequence, the GBI measure- 
ments are only part of the total flux of 3C 273 and have to be 
multiplied by a scaling factor. 

To scale the GBI measurements at 2.25 GHz we use 
the contemporaneous Effelsberg single dish observations at 
2.7 GHz. The nearly flat spectral index of 3C 273 at this fre- 



quency (Tiirler et al. 1999a) allows us to neglect the small dif- 



ference in frequency. We obtain that the GBI fluxes at 2.25 GHz 
have to be multiplied by 1.66 to fit the Effelsberg observa- 
tions. A similar calibration of the earlier GBI measurements 
at 2.7 GHz is not possible due to the lack of contemporaneous 
single dish observations. We therefore let the scaling factor of 



3.1. Typical three-stage evolution 

Following MG85, we consider the synchrotron radiation emit- 
ted just behind a shock front in a cylindric portion of a jet hav- 
ing a radius R across the jet and a length x along the jet axis, 
as illustrated in Fig. |l]. Within this volume, we assume that the 
magnetic field B is uniform in strength and nearly random in 
direction and that the relativistic electrons have a power law en- 
ergy distribution of the form N{E) = K E^", with N{E)dE 
being the number density of the electrons. Measured in the rest 
frame of the quasar, these synchrotron emitting electrons have 
a relativistic bulk velocity (3 = v/c, with c being the speed of 
light, and a corresponding Lorentz factor F = (1 — 
As a consequence, the synchrotron emission observed with an 
angle to the jet axis is Doppler boosted with a bulk Doppler 
factor given hy V = r~^(l — /3 cos6')~^. We further assume 
that the jet opening half-angle is smaller than 6 at any time 
after the onset of the shock and that 9 itself is small enough 
to verify sin (0 + </>)< 1 /F, so that the line-of-sight depth of 
the emitting region is directly proportional to its thickness x 
(Marscher et al. 



1992) 



The observed optically thin flux density and turnover 
frequency of the self-absorbed synchrotron spectrum are 
then given by 



5, oc and (1) 

,2/(s+4) 

(2) 



.^^(s+2)/2 2)(s+2)/2 



' We obtain a best-fit value of 1 . 12 for this factor with the jet model 
presented in Sect. ^. 
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Fig. 1. Sketch of the jet geometry dis- 
cussed in Sect. as observed in the 
rest frame of the quasar. We assume the 
contribution of the inner jet to be neg- 
Hgible with respect to the synchrotron 
emission of the grey region behind the 
shock front. This cyHndrical portion of 
the jet is moving along the jet axis with 
a relativistic speed f3 c 



where K and B are measured in the comoving frame of the 
emitting plasma and S^, v, v^, R, x and V are measured in the 
observer's frame. 

The thickness x of the emitting region behind the shock 
front is the crucial parameter that changes the outburst's evo- 
lution from one stage to the other. If radiative losses are the 
dominant cooling process of the electrons accelerated at the 
shock front, x is approximatively given by 2 v^d icooi, where 
Vroi is the excess velocity of the shock front relative to the 
emitting plasma and tcooi is the typical cooling time of the 
electrons. During the first stage of the shock evolution Comp- 
ton scattering is the dominant cooling process of the elec- 
trons and the Compton limited thickness xi of the emitting 
region is given by xi oc u^^^ B^^^ 2?^/^ iy^^/^, where the en- 
ergy density of the synchrotron photons itph can be expressed 
as Uph cx K (^3^+7 i?*+5)i/8 (Eq. (13) of MG85).0 By com- 
bining these two equations, we obtain^ 



-1/2 



(3) 



The shock evolution enters the synchrotron stage as soon 
as the photon energy density Uph is equal to the magnetic field 
energy density ub — i?^/(87r). During this second stage, syn- 
chrotron radiation is the dominant energy loss process of the 

^ The centre of the emitting region is given by the typical distance 
'^rci tcooi from the shock front reached by the electrons in the upstream 
direction before loosing substantially their energy. The full width x of 
the emitting region is thus about twice that distance. 

^ This calculation assumes that only first order Compton scattering 
is important. The effect of higher order Compton losses would be to 
steep en even more the initial rise of the spectrum (Marscher et al. 



1992) 



■ xi and X2 in Eqs and (Q) are actually expressed in the comov- 
ing frame of the emitting plasma. The transformation of Xi to the ob- 
server's frame would be obtained by dividing these expressions by the 
Lorentz factor T to correct for length contraction and by multiplying 
them by the Doppler factor T) to correct for front-to-back time delay 
(Marscher A.P., private communication). This tricky transformation 
is due to the non-spherical shape of the emitti ng slab , which appears 
rotated in the observer's frame (cf. Marscher 



1987 



Marscher et al. 

1992). However, since the viewing angle 6 is assumed to be smaller 



than about l/F, T> is proportional to T (e.g. Marscher 198C), so that 



electrons. Due to the similarity between Compton and syn- 
chrotron expressions of the cooling time tcooi of the electrons, 
the synchrotron limited thickness X2 of the emitting region is 
simply obtained by replacing Uph in the expression of xi by 
ws, so that X2 becomes 



.T2 0ci3-3/2pl/2^-l/2 



(4) 



During the final stage of the shock evolution synchrotron losses 
become less important than adiabatic expansion losses and the 
shock is assumed to remain self-similar so that 2:3 cx R. 

By substituting x in Eqs (|l|) and ^ by the corresponding 
expression of Xi {i = 1, 2, 3) for each of the three stages i of 
the shock evolution we obtain: 



S^.l oc i?(ll--)/8B(-+l)/8p(-+4)/2^-./2 

S,,2 OC i?2^^(s-2)/2p(s+4)/2^-s/2 

^,,3 OC i?3^5(s+l)/2p(.+3)/2^-(.-l)/2 

OC i?-l/4 5l/4p(.+3)/(.+5) 

1^^^2 OC [i^2^s-lps+3^1/(s+5) 

i.„,3 OC [i?i^B(«+2)/2p(.+2)/2]2/(.+4)^ 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



Finally, by using these six last equations, it is straightforward 
to obtain the expressions for 5m. i = Si,,i{i'ni,i) as 



5,mi oc i?ll/8ijl/8p(3.+10)/(.+5) 
5m,2 OC i?2[^5 52.-5p3.+10jl/(s+5) 
Sm.3 OC [i?2.+13^5 52s+3p3s+7il/(.+4) 



(11) 

(12) 
(13) 



We now assume that each of the quantities K, B and V 
evolves as a power-law with the radius R of the jet and we 
parameterize this evolution as 



K ocR 



-fc 



B ocR 



-b 



V ocR 



-d 



(14) 



If we replace these relations into the expressions for v^^i 
and Sya^i (Eqs (^ to (pj|)), the turnover {v^^, S^) of the self- 
absorbed synchrotron spectrum will also evolve as a power-law 
with R during each stage i of the shock evolution as 



these two corrections cancel out when using proportionalities. 



iyrr,A OC R" 



and Sn 



Rf' 



m.i OC ly^ 



(15) 
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where the exponents rii and fi are given by 

ni = -(6+l)/4 - d (s + 3)/(s + 5) (16) 

712 = -[2fc + 6(.s-l) + d(s + 3)]/(s + 5) (17) 

713 = -[2(fc-l) + (fe + d)(s + 2)]/(s + 4) (18) 
/i = (ll-6)/8-d(3s + 10)/(s + 5) (19) 
/2 = 2- [5fc + fe(2s-5) + d(3.s + 10)]/(s + 5) (20) 
/g = [2s + 13-5fc-fe(2s + 3)-d(3s + 7)]/(s + 4). (21) 



Marscher (1987i shows that even a geometrically thin 



source in a relativistic jet appears inhomogeneous when ob- 
served with a small angle 6 between the jet axis and the line- 
of-sight. The optical depth of the source will therefore depend 
on the frequency, which has the effect of broadening the self- 
absorption turnover and leads to a lower optically thick spec- 
tral index athick than the typical value of +5/2, which holds 
for a homogeneous synchrotron source. To estimate the value 
of the spectral index below the turnover frequency i/,„, we can 
cut the source into many self-similar cylindric portions of the 
jet having a length I <^ x proportional to their radius R and 
chosen small enough for their synchrotron emission to be ho- 
mogeneous. The emitted spectrum of each section will have a 
self-absorption turnover (^'in,S'm) depending on R according 
to Eg. (|l5[) with because in this case / oc R replaces x in 
Eqs ([l]) and (^. The inhomogeneous source behind the shock 
front is therefore expected to have an optically thick spectral 
index of athick = fs /n^ due to the superimposition of the ho- 
mogeneous spectra of the individual sections. However, the fi- 
nite size of the emitting region should limit the frequency range 
over which this flatter spectral index pertains. We therefore ex- 
pect a spectral break at a frequency i^h (< i^m), at which the 
spectral index athick = fs/n^ returns to its homogeneous value 
of +5/2. 

A high-frequency spectral break is also expected in the op- 
tically thin part of the spectrum due to a change in the elec- 
tron energy distribution induced by synchrotron and/or Comp- 
ton losses. In the case of continuous injection or re-acceleration 
of electrons suffering radiative losses, the optically thin spectral 
index athin is expected to steepen by a value of —1/2 above a 
frequency i^b (Kardashev 1962). The break frequency i^b of the 



spectrum in the observer's frame is related to the break energy 
iJb of the electron energy distribution as i-'b oc BVE^ (e.g. 



Marscher 1980). For our analysis at submillimetre-to-radio fre- 



quencies, the evolution of i/b with the jet radius R is only rele- 
vant during the final stage of the shock evolution. For adiabatic 
expansion in two dimensions, the energy E of th e electrons de- 
creases with i? as i? oc R^^^^ (e.g. Gear 1988). By using this 
relation and Eq. ([l4|), we find that the evolution of the break 
frequency I'h during the adiabatic stage is given by 



(22) 



z/b,3 OC 7?"*= with nyy = -{A/3 + b+d) . 

Until now, we related the spectral turnover {i/^, Sai) and 
the high-frequency spectral break rib to the radius R of the 
emitting region. But since we are interested in the temporal 
evolution of these quantities, we have to relate R to the ob- 
served time t after the onset of the outburst. According to the 



geometry of Fig. |l]and the basic principles of superluminal mo- 
tion (e.g. Pearson & Zensus 1987), we obtain 



(l + z)sin^ 

Mapp ^ 



(3cTV 



(23) 



where /3app is the apparent transverse velocity of the source 
in units of c and L measures the distance along the jet axis in 
the rest frame of the quasar. To continue working only with 
proportionalities, we are forced to assume that Lon is small 
with respect to L during the major part of the shock evolution, 
so that L — Lon tends rapidly towards L during, or just after, 
the Compton stage. Under this assumption and by considering 
that (3 r V is pr oportional to for T > 1 and 9 < l/T (e.g. 
Marscher 1980), we have t oc L. Finally, by parameter- 



izing the opening radius R of the jet with the distance L as 
R oc L'' and by remembering that V oc R^'^ (Eq. ([l4l)), we 
obtain 



t(xV-^ R^^'' oi RP with {2rd+l)/r. 



(24) 



With this relation, the three-stage evolution of Eq. ( |15| ) can now 
be expressed as a function of time t rather than radius R, as 



t'ma oc t^* with Pi = ni/p 
Sn,.i oc with -fi = fi/p 



2i 



(25) 



where (3i and 7^ are defined as in Paper 1 and depend now on 
the five exponents s, r, k, b and d via the Eqs ( [T^ to (^T]) for 
and fi (i = l,2,3). 

3.2. Specificity of individual outbursts 

We described in the previous section how the spectral turnover 
{i^m, Sm) of the emitted synchrotron spectrum behaves accord- 
ing to the physical properties of the jet and during each of the 
three stages of the shock evolution. What we describe here is 
how the transitions from one stage to the other are influenced 
by the physical quantities at the onset of the shock. The first 
transition from the Compton to the synchrotron stage is char- 
acterized by the condition Uph = ub, which corresponds to 
xi — X2 and can be expressed by the proportionality 



R 



-(s+5) 
1|2 



'1|2 -°1|2 



3(s-3) 



(26) 



where /ii|2 = i^(ti|2), Ri\2 = R{ti\2) and 3^2 = -B(ti|2), 
with the subscript 1|2 referring to the transition from the first 
to the second stage of the shock evolution at a time ti|2 after 
the onset of the shock. 

A similar expression can be derived for the second transi- 
tion from the synchrotron to the adiabatic stage by imposing 

3/9 1/2 1/2 

that X2 — X3, which corresponds to i?2|3 oc B^^^' "2^213 ^^213 ' 
where j/2|3 = i^m(i2|3)i with the subscript 2|3 referring to the 
transition from the second to the third stage of the shock evo- 
lution. By replacing 1^2 1 3 in the condition above by the propor- 
tionality of Eq. (^ for fm, 2(^213) - or equivalently by that of 
Eq. ^) for i^m, 3(^213) - we obtain 

-"213 « ■"-2|3-fc'2|3 •'^213 • 



(27) 
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We can now further parameterize the proportionalities of 
Eq. (|l4|), as 



K cx Kon R 



-k 



B (xBonR 



-b 



V OC Von R 



-d 



(28) 



where the subscript "on" refers to the onset of the shock. By 
including these relations expressed at the particular radii R — 
i?i|2 and i? = i?2|3 into Eq. ( |2^ and Eq. (p7|), respectively, we 
obtain 



r[]^^ oc Bfi'-^^ and i?! 



on on 



2|3 ^ ^on ^on ^on ' (29) 



where the exponents ^112 ™d ^213 are given by 



Cl|2 
C2I3 



8fc + 36(s-3)- (s + 5) 
fc + 6(2s + 7)-d-(s + 5). 



(30) 
(31) 



The expressions of Eq. ( |29| ) relate the radii i?i|2 and i?2|3 
of the jet at which the two transitions 1|2 and 2|3 occur to the 
values of the physical quantities K, B and V at, or just after, 
the onset of the shock. Any change of these quantities at the 
onset of the shock will therefore have the effect of displacing 
the position in the jet (via L cx R^^^) at which the two tran- 
sitions occur A stronger magnetic field Bon, for instance, will 
have the effect to prolong the synchrotron stage, by making it 
start slightly further upstream in the jet and finish much further 
downstream for a typical value of s w 2 and comparable values 

of Ci|2 andC2|3- 

The values of Kon, Bon and Von, will also influence the 
times ii|2 and t2|3 after the onset of the shock at which the 
two transitions are observed, as well as the frequencies i'n2 
and i^2\3 and flux densities 5i|2 and 5*213 of the self-absorption 
turnover (i/m, 5111) at these times. By including the proportion- 
aUties of Eq. into Eq. ( ^4| ) for the observed time t and into 
Eqs and (|l2|) for the evolution during the synchrotron stage 
of turnover frequency I'm, 2 and the corresponding flux density 
Sm.2, we obtain for the first transition 



ii 19 cx; Von Rl\2 



^l|2 
J^l|2 OC 



(pit' Kl^Bl-^^'^'' R-^ 



1|2 



5i|2 OC {Vli+''Kl^B^o!.^'Y'^'^'^R{l 



(32) 
(33) 
(34) 



Exactly the same expressions apply also to the transition 213, 
because both transitions are related to the synchrotron stage. ^ 
If we now substitute the expressions of Eq. (29) for i?i|2 
and i?2|3, we obtain the proportionalities which define the place 
in the {S, v, i)-space (cf. Fig. ||) where the transitions from one 
stage to the other occur according to the values of Kon, Bon and 
Von- These proportionalities can be expressed as logarithmic 
shifts A Ig P = Ig P — (Ig P) from an average value (Ig P) of 
the parameter P, which stands either for Kon, Bon and Von or 



We would have obtained other proportionalities by using the ex- 
pressions of the Compton stage for the first transition or those of 
the adiabatic stage for the second transition. However, the end re- 
sult of Eqs ( p6[ ) and ( |37| ) would have been the same, because 1/^2 = 
J^m,i(ti|2) = ;^m,2(ti|2). and similarly for Si|2, z^2|3 and 5'2j3- 



Table 1. Expressions of the parameters U, V and W of Eqs ( 
to and of Eq. for the logarithmic shifts A Ig P of the 
parameters P in the first column. These expressions character- 
ize the effect of varying the physical quantities K, B and V at 
the onset of the shock. 





U 


V 


w 


ii|2 


8p 

Cl|2 


3(3-3)p 

Cl|2 


-2 


^213 


P 

C2I3 


(23+7)p 
C213 


'>2|3 


1^112 


2 1 8n2 
s + 5 ^ Cl|2 


s — 1 1 3(s — 3)n2 
3 + 5 Cl|2 


3 + 3 
3 + 5 


1^213 


2 1 112 
3+5 ^ C213 


s-1 1 (2s + 7)n2 
3 + 5 C2|3 


3+3 _ n2 


3+5 C213 


'S'l|2 


5 1 8/2 
3 + 5 ^ Cl|2 


2s-5 1 3(s-3)/2 
3 + 5 Cl|2 


3s+10 

3 + 5 


'S'213 


5 1 /2 
3 + 5 ^ C2|3 


2s-5 1 (23+7)/2 
3 + 5 C2|3 


33+10 /2 
3+5 C2I3 


'^b,2|3 


"b 
C2I3 


I 1 (23 + 7)nb 
C2I3 


1 _ "b 
C2I3 



for S, V and t. The set of equations relating these shifts for the 
first transition can be written as 



-Vt„2Algi?o 



K„2AlgSo 



Wt,,,/^\gVon (35) 
W.^I^AlgPon (36) 
-VFs,,2AlgPon.(37) 



Algti|2 = C/t,|2Algi^, 
Algj.i|2 = C/.,|2Algi^, 

Alg5i|2 = [/5i|2 Algi^on + l^Si|2 AlgP, 

Similar expressions can be written for the second transition by 
replacing 1|2 in Eqs (^sj) to ( ^7| ) by 2|3. The 2x9 expressions 
of the parameters U, V and W obtained for both transitions are 
displayed in Table ^ There is a great symmetry among these 
parameters, which reflects the fact that the two equations 

Algi/2|3 - Alg:/i|2 = /?2(Algi2|3 - Algii|2) and (38) 

Alg52|3-Alg^l|2 = 72(Algi2|3- Algti|2) (39) 

must be verified, in accordance with Eq. (p5[). 

Finally, we can also derive the influence of Kon, Bon and 
Von on the high-frequency spectral break ^'b,2|3 = ^'b(i2|3), 
which is given by ^'b,2|3 oc Von Bon -R213 (''^- '^l- (p2|)) and 
corresponds to the logarithmic shift 

Alg l^b,2|3 = C^^b Algi^on + Kb AlgPo, 

with the expressions of Table |l] for 
4. Parameterization 

In the previous section, we described in detail the physical ba- 
sis of the light-curve decomposition that we present in Sect. |[ 
Here, we are more concerned by the practical aspect of the de- 
composition and particularly by the proper definition of the 85 
parameters used by the fit. 

4.1. The synchrotron spectrum and its average evolution 



-W^.bAlgPo 



andM^^j,. 



(40) 



According to Sect. 3.1, the evolution of a shock wave in a jet 
follows three distinct stages to which correspond different evo- 
lutions of the self-absorption turnover of the synchrotron spec- 
trum emitted by the plasma behind the shock front. The evolu- 
tion of the turnover {vni, Sm) during each of the three stages is 
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fully determined through Eq. ( p5| ) by only five indices, which 
are the three indices fc, b and d defined by Eq. (|l4|), the index 
s of the electron energy distribution N{E) and the index r of 
the relation R oc L"^ defining how fast the jet opens with dis- 
tance L. In the model presented in Sect. ^, we impose the value 
of d to be zero and leave the four other indices as free param- 
eters of the fit. They are constrained by the slopes /3i and 7^ 
(i = 1, 2, 3) of the outburst's evolution in the (Ig S, Ig h', \gt)- 
space (cf. Fig. ^ and are the most interesting parameters of 
the light-curve decomposition, because they describe physi- 
cal properties of the jet. The two points (ti|2, i^i|2i '5'i|2) and 
(i2|3i 1^213, 5'2|3) in this space at which the transitions 1|2 and 
2|3 from one stage to the other occur are defined by only four 
of these six quantities, because the slopes (32 and 72 can be used 
to determine 1^2 1 3 and iS'213 from the four other values. 

Until now, we have used eight parameters, namely s, r, k, 
b, lg<i|2, lgi2|3^ lgf^i|2 and lgiS'i|2, to fully characterize the 
path in the {S, u, t)-space followed in average by the maxi- 
mum (j/,n, S'm) of the self-absorbed synchrotron spectrum. The 
shape of this spectrum is defined by the general expression (cf. 
Paper I) 



where Tn 



l-exp(-Ti„ (zy/zy^ "thick) 

1 - exp (-Tm) 



(41) 



1 - 



- 1 



g^^^ ^ i , is a good approximation of 
the optical depth Ti, at the turnover frequency and Si^ is the 
real maximum of the spectrum.^ We choose to impose the value 
of the optically thick spectral index athick to /a / nj, according 
to the discussion in Sect. 3.1 . The optically thin spectral index 
Qfthin is also determined by the model and has the value — s/2 
during the first two stages before flattening to — (s — 1)/2 at 
the transition 2|3 to the adiabatic stage, according to Eqs (^ to 
(^. This flattening by Aathin = +1/2 is assumed to end at 
the time t2|3 of the transition, but since it cannot be instanta- 
neous we are forced to use one parameter to define the time tf 
when it starts. Because this spectral change is expected to begin 
slowly before accelerating until the transition 2|3, we choose to 
describe it with a logarithmic expression of the time t, as 



athin(*) = 



1 igftAf) 

2 lgfe|3Af) 



for ti <t< t2\3 ■ 



(42) 



According to the considerations of Sect. 3.1 
low the spectrum defined by Eq. (STl) to have a high 



we further al- 
and a 



low-frequency spectral break as shown in Fig. ^. Since the ex 



act shape of the breaks is difficult to assess (cf. Marscher 1977 



Band & Grindlay 1985), we choose to keep them sharp, rather 



than smooth them arbitrarily, because sharp breaks have the ad- 
vantage to define precisely the frequency at which they occur 
The spectrum therefore simplifies to 6*11 (z^/z^h)^'^ at the lowest 




Frequency 

Fig. 2. The shape of the synchrotron spectrum (solid line) as- 
sumed to be emitted behind the shock front during the final 
stage of the shock evolution. We show the effect of adding two 
spectral breaks to the simpler spectrum (dashed line) used in 
Paper I 



frequencies {ly < i^i^) and to 5b (i^/t'b) at the highest fre- 
quencies {1/ > Vh), where Sh = Si,{h'h) and Sh = S^{i^h) are 
calculated by using Eq. (|4l|). The frequency ratio t/h/ I'm of the 
low-frequency break to the spectral turnover is a free parame- 
ter of the fit and is assumed to remain constant throughout the 
outburst's evolution. This ratio might actually slightly increase 
with time due to the increase of the thickness x of the emitting 
region behind the shock front, but, for simplicity and because 
the effect might change from one stage to the other, we do not 
take this possible increase into account. 

Since we do not extend the light-curve decomposition up 
to infrared frequencies, we do not need to consider a high- 
frequency spectral break during the two first stages of the 
outburst's evolution, which would have the effect to further 
steepen the optically thin spectral index atinn to — (s + l)/2. 
For frequencies above the break frequency i^b, athin remains 
therefore always at its value of — s/2, even when it is flat- 
tening at lower frequencies. The evolution with time of the 
break frequency z^b is given by Vh — t'b,2|3 (V^2|3)"''^''^ where 
z^b.2|3 = '^h{t2\3) is a free parameter of the fit and the ratio 
nb//0 is fixed by the values of r, b and d through Eqs ( p2| ) and 



Note that 5*^ as defined in Sect. 3.1 



the extrapolation S^*^'" = 5'J'^'"(i-'m) down to i/m of the optically 
thin spectrum. The real maximum Sm is related to 5^*^'" by 5m = 
Sm^'°(l — exp (— rm))/''m, but this distinction is not important for 
this work, which is based on proportionalities. 



is not the real maximum, but 4.2. Characteristics of different outbursts 



The 11 parameters defined above suffice to describe fully the 
spectral evolution of the average outburst shown in Fig. |} The 
behaviour of the spectral turnover during each of the three 
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stages of the evolution will be exactly the same for all individ- 
ual outbursts. What can change from one outburst to the other 
is the position in the (Ig S, Ig ly, Ig i)-space where the two stage 
transitions occur In Paper I, the two transition points were not 
allowed to move separately and therefore all outbursts were 
self-similar in the sense that they all had exactly the same evo- 
lution pattern in the (Ig 5, Ig i^, lgt)-space. Here, the distance 
between the two transition points can change and hence make 
the synchrotron stage relatively longer or shorter. 

Although this new effect makes the light-curve modelling 
more complex, we do not introduce a further free parameter per 
outburst to describe it. Instead, we now give a physical origin 
to the three logarithmic shifts A Ig S, Algiy and A Ig < in flux 
density S, frequency i/ and time t of Paper I, by replacing them 
by shifts of the values at the onset of the shock of the factor 
K in the expression of the electron energy distribution N{E), 
the magnetic field B and the Doppler factor V. The change of 
the three physical quantities Kon, ^on and X'on results in dif- 
ferent shifts in flux, frequency and time for the two transitions 
1|2 and 2|3 according to Eqs ( ^5[ ) to (p7[). The shifts used in 
Paper I appear therefore to be quite unrealistic, because any 
such displacement in the (Ig S, Ig Ig i)-space will also result 
in making the synchrotron stage longer or shorter. 

Actually, the physical link between the shift of an outburst 
and its shape is able to solve some problems we had in Pa- 
per I (cf. Sect. H), but does not improve the overall quality 
of the light-curve decomposition. To obtain a satisfactory fit, 
we have to add two outbursts, one in 1988 and one in 1991. 
Three other outbursts were added to extend the modelling of 
the light-curves to the most recent observations included here 
(cf. Sect. We end up with seventeen outbursts, each charac- 
terized by four parameters, namely the date of its onset Ton and 
the three logarithmic changes A Ig i^Ton, A Ig Son and A Ig 2?on 
of the physical quantities K, B and V at the onset of the shock. 
We therefore need a total of 17 x 4 = 68 parameters to model 
the specificity of the seventeen individual outbursts. 

4.3. Contributions from the outer jet and previous outbursts 

The outer jet and especially its terminal hot spot called 
3C 273A contributes significantly to the observed flux density 
at low radio frequencies, but this contribution can be consid- 
ered as constant during the 20 years of this work. Th e spec trum 
of the hot spot 3C 273 A is shown by Conway et al. ( 1993 ) and 
corresponds roughly to a power-law with a spectral index a of 
0.85 above ^ 1 GHz, as described in Tiirler et al. ( 1999a ). As 
in Paper I, the fit of the light-curves by the series of outbursts 
is done above this constant power-law contribution. 

Another important contribution at low frequencies comes 
from the superimposed decays of the outbursts peaking before 
1979. In Paper I, we modelled this contribution by an exponen- 
tial flux decay, which was appropriate because the decline of 
the individual outbursts was also modelled by an exponential 
decay in the light-curve approach, so that the characteristic e- 
folding time was well constrained at all frequencies. It would 
be unsatisfactory to take here the same exponential decays, be- 



cause the declining phase of a three-stage evolution implies that 
the final light-curve decay of each outburst is a power-law and 
not an exponential, as illustrated by the linear ending of the 
logarithmic light-curves of Fig. 

The observed flux decay Sl{t) = S^{h'*,t) with time 
i at a given frequency i^* can be written as S*{t) — 
<S'm(t)(j^*/fm(^))"""", where the equahty holds strictly if 
is actually the extrapolation S^^^™ of the optically thin spec- 
trum (cf. footnote ||). During the final stage of the outburst's 
evolution, the time dependence of this quantity can be ex- 
pressed as 5m(t) = 5'*,(i/im)''' (cf- Eq. (H)), where = 
tm ( ) is the time after the onset of the outburst when the spec- 
tral turnover {ly^, S^) passes at the frequency i/*, i.e. the time 
t for which i/ni(0 = '^*- The corresponding equation for the 
turnover frequency is t',n(^) = '^m(V^m)''^- With these two 
equations, the light-curve decay at a given frequency ly* during 
the adiabatic stage is given by 

/ J. \ 73-/33 Qthin / * \ l/ft 

5:(t)ocy with 4-^.,3(— j .(43) 

The quantity ij^ acts as a characteristic timescale of the power- 
law decay at a frequency i/* and its frequency dependence de- 
scribes how this timescale increases towards lower frequencies. 
In Eq. (p3l), athin can be replaced by — (s — 1)/2, if we do not 
consider a possible effect of the optically thin spectral break, 
which arises only at later times, once the break frequency i^b 
has reached the frequency i^* . 

On this basis, we can model the contribution at a given fre- 
quency ly of the superimposed decays of the outbursts peaking 
before 1979 by the decay of an hypothetical outburst having an 
amplitude Ao{iy) at the date Tq = Ton + tra{v) = 1979.0 and 
then decaying as 

/ T^_iq7q\ 73+/33(^-i)/2 
S,{T)^Ao{v)[l + —^] , (44) 

where we have used Eq. ( p3| ) with t replaced by T — Ton = 
tm{i^) + T — 1979 and we have introduced a factor /i (> 1), 
which is frequency independent. There is no physical justifica- 
tion for introducing this factor /i, but it was found to be a sim- 
ple way to take into account the fact that the superimposition 
of several outbursts with different onset dates does decay with 
a longer characteristic timescale as a single outburst. For the 
whole range of frequencies covered by the thirteen light-curves 
considered in this work, Ao(i^) in Eq. (Q) is the millimetre-to- 
radio spectrum of 3C 273 in 1979.0. This spectrum is modelled 
by a cubic spline parameterized at the four frequencies defined 
by Ig (i^/GHz) = 0.3, 0.5, 1.0 and 1.5 and is extrapolated to 
higher frequencies. We choose this low frequency range, be- 
cause the influence of the spectrum Ao{i') is the greatest in the 
radio domain and becomes negligible at submillimetre wave- 
lengths (cf. Fig. H). The value and the frequency dependence of 
the timescale t^{i') is completely determined by the right part 
of Eq. ( p3| ) with v* replaced by v, so that the contribution from 
the outbursts prior to 1979 is only modelled by five free param- 
eters, namely the values of Aq{i') at the four above-mentioned 
frequencies and the frequency independent parameter ^. 
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Fig. 3. Ten out of thirteen submillimetre-to-radio light-curves of 3C 273 (grey points) decomposed into a series of seventeen 
synchrotron outbursts. The best-fit hght-curve (solid line) is the sum of the constant emission from the outer jet (long-dashed 
line), the global decay of all outbursts peakingbefore 1979 (short-dashed line) and the seventeen superimposed outbursts (dotted 
lines) starting at the dates Ton given in Table |3] 



A final parameter of the fit is the GBI scaling factor men- 
tioned in Sect. ^ and therefore the total number of free parame- 
ters used to adjust the 5234 observational points in the thirteen 
light-curves is 11 + 68 + 6 = 85. We fit all light curves simulta- 
neously, but usually with not more than ten of the 85 parameters 
at a time. The four indices s, r, k and b are adjusted together, 
followed by the seven other parameters defining the spectral 
evolution of the average outburst. The specificity of individ- 
ual outbursts is adjusted by a series of fits concerning only two 
or three outbursts at a time. Finally, we fit together the six re- 
maining parameters, which define the initial flux decay and the 
GBI scaling factor. This procedure is repeated many times un- 
til the values of the parameters converge to their best-fit values. 
In total, several hundreds iterative fits are needed to obtain the 
hght-curve decomposition of Fig. ||. 

5. Results 

We present here the results obtained by fitting the 
submillimetre-to-radio light-curves of 3C 273 according 
to the jet model described in Sect. ^ and with the set of param- 



eters defined in Sect. ^ We choose to focus in this section on 
the results obtained by assuming that the synchrotron emitting 
plasma behind the shock front moves with constant speed 
and direction, so that the Doppler factor V remains constant 
during the outburst's evolution. The index d characterizing the 
decrease of T> with jet radius R (cf. Eq. (|l4|)) was therefore 
fixed at a value of zero. All following figures and tables 
correspond to this particular jet model. In Sect. ^+| , we will 
however also discuss the results of an alternative jet model, in 
which the Doppler factor is free to vary, but we impose the jet 
to be conical. 

The best-fit decomposition we could achieve with the jet 
model having a constant Doppler factor V is illustrated in 
Fig. ^ We note that the light-curves which are the best repro- 
duced by the model are the radio light-curves, and especially 
those at 15 GHz and at 8 GHz. In the 22 GHz hght-curve, flie 
rapid variations from 1981.0 to 1983.5 are not well described 
by the model, and at 37 GHz, there is mainly a problem with the 
prominent peak of 1 99 1 .5 . At 90 GHz, it is striking to note that 
the model reproduces quite well the features of the observed 
light-curve, but cannot reach a sufficiently high flux level. The 
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Table 2. Best-fit values of the parameters defined in Sect. ^ and 
corresponding to the evolution of the typical outburst shown in 
Fig. H Other related quantities defined in Sect. |3] are also dis- 
played here and are distinguished from the adjusted parameters 
by the symbol f 



Par. 


Val. 


Par. 


Val. 


Par. 


Val. 


Par. 


Val. 


s 


2.05 


r 


0.82 


k 


3.03 


h 


1.58 


*1|2 


0.07 yr 




442 GHz 




14.2 Jy 




12.7 


h\3 


1.06 yr 




36.7 GHz 


"^213 


16.0 Jy 


C2II 


13.5 


ti 


0.75 yr 




16.3 THz 




0.40 


M 


4.13 




-0.64 




-1.09 




-1.73 


-J 


-2.91 




+ 1.18 




+0.05 




-1.53 




+ 1.22 




-0.53 




-0.89 




-1.41 






7^ 


+0.96 


7j 


+0.04 




-1.25 






71 //3/ 


-1.83 




-0.05 




+0.88 
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Fig. 4. Model light-curves of the average outburst in 3C 273 
at different frequencies spaced by 0.1 dex (grey lines). The 
six highlighted light-curves are at frequencies i' defined 
by lg(i//GHz) — 3.0, 2.5,..., 0.5, in order of increasing 
timescales 



same problem is also apparent at a wavelength of 2 mm and it is 
only at 1 . 1 mm that the model light-curve is again on average at 
the same flux level than the observations. Finally, at the highest 
frequencies, the model light-curve does not decrease enough 
between the outbursts to adjust the lowest flux observations. 

This general discrepancy at millimetre and submillimetre 
frequencies is probably responsible for the relatively high re- 
duced value of Xrcd = X^/d-O.f. = 20.5 obtained for the 
overall fit as compared to the values quoted in Paper I. Apart 
from this global problem, which will be discussed in Sect, ^j] , 
we note that the detailed structure of the light-curves is now 
better reproduced by the set of outbursts than it was in Pa- 
per I. In particular, the fast decay around 1985 is now more 
satisfactorily accounted for by the model. The rapid variations 
of the flare peaking near 1996.0 are also better described by 
the new parameterization. If we compare the high-frequency 
light-curves of Fig. ^ with those shown in Fig. 7b of Paper I, 
we see that the improved description of this outburst starting in 
1995.8 is achieved with a faster decline just after the peak of the 
high-frequency light-curves. This faster decline results from 
the introduction of the new parameter <t, which defines the time 
when the optically thin spectral index begins to flatten and from 
the fact that this time is delayed with respect to the time <i|2 of 
the first stage transition. This flattening, occurring now just at 
the end of the synchrotron stage, is in very good agreement 
with the shock model of MG85. The above-mentioned better 
description of the decay around 1985 is not due to the same 
effect. A satisfactory fit of 3C 273's behaviour at this epoch 
could indeed only be achieved once we took into account the 
possibility that the stage transition points 1|2 and 2|3 could be 
shifted separately as it is the case here (cf. Sect. 4.2 ). The ap- 
propriate shape of the last three outbursts peaking before 1985 
was obtained by reducing very much their intermediate syn- 
chrotron stage, as shown by the values in the last column of 
Table H 

The best-fit values of the adjusted parameters and the cor- 
responding values of other interesting quantities are displayed 



in Table ^ These values define the three-dimensional evolution 
of the average outburst in 3C 273, as shown in Fig. |] We note 
that the decline of the turnover flux 5^ with turnover frequency 
z^ni is now flatter than obtained with the three-stage approach 
of Paper I. Instead of "fz/ = 1.14, we now find 737/33 = 0.88. 
This difference could be due to the introduction of the high- 
frequency break in the spectrum (cf. Fig. Another difference 
is the steeper rise of with by nearly a factor 2 (71 / (3i 
is now of —1.83 instead of —0.99). This might be due to the 
introduction of the low-frequency spectral break, but it must 
be noticed that this initial part of the evolution is the less con- 
strained stage of the outburst's evolution. Its behaviour could 
therefore be imposed by the best-fit values of s, r, k and b, 
which are more strongly constrained by the two other stages of 
the evolution. Nevertheless, both the steeper rise and the flatter 
decline obtained here are in better agreement with the original 
outburst's evolution proposed by MG85. 

The value of 2.05 that we obtain for the index s of the elec- 
tron energy distribution N{E) corresponds to an optically thin 
spectral index of atinn = —{s — l)/2 w —0.53 for the final 
adiabatic stage of the outburst's evolution. These values for s 
and athin are slightly higher than those obtained in Paper I, but 
are still well below the values of s ~ 2.4 and atinn = — 0.7 
measured by MG85 for the strong flare of 1983. The parameter 
K of the electron energy distribution N{E) is found to evolve 
with the width of the jet R approximately as if oc which 
is just a bit steeper than the decrease K cx i?~2(s+2)/3 ^ 



expected if the jet flow was adiabatic (e.g. Gear 1988). The de 



crease of the magnetic field B with R is found to be between 
the two extreme cases of i? oc R^^ expected if B was paral- 
lel to the jet axis and B oc R~^ expected if B was transverse 
to the jet (e.g. Begelman et al 
b — 1.58 suggests that the par 
B± components of the magnetic field are roughly equal. Fi- 
nally, the best-fit value of r = 0.82 obtained with this model 
having a constant Doppler factor V suggests that the inner jet 
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and the perpendicular 
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Fig. 5. a-d. Logarithmic evolution of 
the average outburst in 3C 273. Panel d 
shows the three-dimensional represen- 
tation in the (Ig S, Ig Ig i)-space. The 
other panels show the three Cartesian 
projections of this surface: a logarith- 
mic light-curves at different frequen- 
cies spaced by 0.2 dex; b synchrotron 
spectra at different times spaced by 
0.2 dex; c contour plot in the frequency 
versus time plane. The thick solid line 
is the path followed by the turnover 
of the spectrum. This line is a bro- 
ken power-law having three different 
indices in each panel. From left to right, 
these indices are 71, 72, 73 in panel a; 
Is/Ps, 72//32, 7i//3i in panel b; and 
Pi, P2, P3 in panel c, with the val- 
ues given in Table |[ The dashed line 
is the line connecting the peak fluxes 
of the light-curves at different frequen- 
cies. The points show for each of the 
seventeen outbursts the position that 
would have the two stage transitions 
1|2 (triangles) and 2|3 (circles) as a 
consequence of the logarithmic shifts 
A Ig Kon, A Ig Bon and A Ig 2?on given 
in Table I 



of 3C 273 is not really conical, but tends to open slightly less 
with distance L along the jet axis (cf. Fig. 

The individual characteristics of the outbursts are pre- 
sented in Table |[ They are expressed by the logarithmic shifts 
A Ig P = Ig P — (Ig P) from an unknown average value (Ig P) 
of the physical quantities Kon, Bon and Von at the onset of 
the shock. The dispersions a of the seventeen values in the 
columns of Table || are cr w 0.40 for A Ig Kon, o ~ 0.29 for 
A Ig Bon, and a « 0.09 for A Ig Von- We note that the Doppler 
factor V is the quantity with the smallest relative changes from 
one outburst to the other. To study the possibility that the log- 
arithmic shifts of the three physical quantities Kon, Bon and 
T>on are correlat ed, w e applied a Spearman rank-order test 
(e.g. Press et al. 1992) to these three data sets. We find out 



that there is a very significant anti-coiTelation between the val- 
ues of AlgA'on and those of AlgI?on- The probability that 
a stronger anti-correlation could occur by chance is less than 
6 10^^. A weaker trend, with a probability of non-correlation 
of less than 2 %, suggests that A Ig T>on is positively correlated 
with A Ig Bon- On the other hand, no significant correlation is 
found between A Ig Kon and A Ig Bon- These relations are dis- 



cussed in Sect. 5.2 



6. Discussion 

In Paper I, we discussed the relationship between the onset of 
an outburst and the ejection of a new VLBl component in the 
jet and we concluded that the identified outbursts do fairly well 
correspond to the VLBl knots. The model presented here does 
not change this conclusion and therefore we postpone a further 
discussion of this relationship until more recent VL Bl m easure- 
ments of 3C 273 are published. We focus, in Sect. 6A_ , on the 
global properties of the inner jet, before discussing, in Sect. |6.2| , 
the peculiarities of individual outbursts. 



6.1. On the global properties of the inner jet 

We have shown in Fig. || the decomposition of a series of light- 
curves. It is, of course, also possible to show the corresponding 
decomposition of the observed spectrum at different epochs. 
As an example, we give in figure Fig. ^ the spectral dissec- 
tion of the two most different submillimetre-to-radio spectra 
observed in 3C 273 during the last 20 years. Until now, such 
a spectral dissection could only be achieved with VLBl obser- 
vations and the best example of this is probably the decompo- 



sition by Marscher (1988) of the radio spectrum of the rela- 



tively high redshift quasar NRAO 140. Actually, the great sim- 
ilarity between the dissection of NRAO 140, based on a sin- 
gle epoch VLBl image, and the dissection of Fig. ^, based on 
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Fig. 6. a and b. Spectral dissection of 
3C 273's submillimetre-to-radio emis- 
sion in the 2-2000 GHz range at the 
epochs of highest (a) and lowest 
(b) flux levels in the submillimetre 
range during the last 20 years. These 
epochs coiTespond to the famous spec- 
tra of 3C 273 pu blished by Robson 
et al. ( |l983i |l986| ). The model spec- 
trum (solid line) is the sum of the 
constant outer jet contribution (long- 
dashed line), the slowly decaying con- 
tribution from the outbursts peaking be- 
fore 1979 (short-dashed line) and the 
following individual synchrotron out- 
bursts (dotted lines) evolving from high 
to low frequencies 



Table 3. Individual characteristics of the seventeen outbursts 
starting at the dates Ton given in the first column. The three 
next columns display, respectively, the logarithmic shifts of the 
normalization Kon of the electron energy distribution N{E), 
of the magnetic field Bon and the Doppler factor Pon, with 
all three quantities evaluated at the onset of the shock. The two 
last columns show the effect of these shifts on the relative max- 
imum flux reached by each outburst, as expressed by Alg S'213, 
and the relative duration of the synchrotron stage, as expressed 
by lg(t2|3/ii|2)- In Fig- ^a, these two last quantities character- 
ize, respectively, the flux density level of the circles and the 
distance between them and the corresponding triangles 



Ton 


AlgJ^on 


Alg Bon 


AlgOon 


Alg52|, 


lg(i2|3Al|2) 


1979.56 


+0.08 


-0.21 


-0.07 


-0.09 


0.90 


1980.41 


+0.42 


+0.22 


-0.02 


+0.23 


1.21 


1982.30 


+0.49 


-0.58 


-0.09 


+0.19 


0.14 


1983.05 


+0.67 


-0.43 


-0.07 


+0.36 


0.21 


1984.02 


+0.57 


-0.41 


-0.07 


+0.27 


0.30 


1986.32 


-0.18 


+0.08 


+0.06 


-0.01 


1.43 


1987.90 


+0.17 


+0.19 


-0.05 


+0.00 


1.34 


1988.19 


-0.48 


+0.05 


+0.06 


-0.20 


1.60 


1990.34 


-0.04 


+0.56 


+0.07 


+0.10 


1.95 


1990.92 


-0.24 


+0.19 


+0.09 


+0.01 


1.61 


1991.12 


+0.28 


-0.32 


-0.21 


-0.26 


0.63 


1993.33 


-0.39 


+0.18 


+0.09 


-0.10 


1.69 


1994.64 


-0.50 


-0.01 


+0.04 


-0.26 


1.53 


1995.84 


-0.15 


+0.09 


+0.08 


+0.08 


1.42 


1996.74 


+0.11 


+0.21 


-0.03 


+0.00 


1.40 


1997.35 


-0.11 


+0.08 


-0.00 


-0.09 


1.39 


1998.30 


-0.71 


+0.10 


+0.13 


-0.22 


1.81 



multi-wavelength variability observations, gives strong support 
to the idea that VLBl knots, observed in the jet, are physically 
linked to the outbursts seen in the light-curves. 



The kind of spectral shape shown in Fig. with a rel- 
ative maximum of the spectrum in the milli metre range, was 
observed in several blazars by Brown et al. ( 1989 ) and inter- 
preted as the evidence of two synchrotron components con- 
tributing to the spectrum. The two decompositions of Fig. ^ 
suggest that these two components are not of distinct nature, 
but that the flaring component progressively becomes part of 
the more slowly varying component as it evolves towards lower 
frequencies. It means that, apart from the contribution of the 
outer jet, the quiescent radio emission of 3C 273 can be entirely 
attributed to the superimposed decays of earlier outbursts. 

We have assumed in this work that the emission of the in- 
ner jet is negligible outside of a small region containing the 
shocked plasma just behind the shock front (cf. Fig. |l|). It 
means that we have implicitly made the hypothesis that the 
emission of the unresolved VLBl core is due to the emission 
of the one or two most recent outbursts observed during the 
earlier phases of their evolution. This is a rather extreme as- 
sumption, since it is usually believed that the base of the jet 
is also contributing to the emission, because of the higher elec- 
tron and magnetic field density at the narrow end of the jet. The 
present work shows that the main features of submillimetre-to- 
radio emission of 3C 273 can be understood as being only due 
to Shockwaves. There are, however, some indications that there 
is still some place left for an underlying contribution of the in- 
ner jet. The presence of an inner jet contribution is suggested 
by the relatively poor fit at millimetre and submillimetre wave- 
lengths and especially by the fact that it is not the shape of the 
variations which is the greatest problem, but much more the av- 
erage flux level reached by the model, which is too low at some 
frequencies and too high at others (cf. Sect. ||). The idea behind 
this is that the model seems to make a compromise between the 
opposite requests of, on the one hand, the millimetre observa- 
tions asking to flatten the optically thin spectral index to have 
more flux during the slow decline of the outbursts and, on the 
other hand, the submillimetre observations asking to steepen 
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this index, in order to decrease the flux level at the end of the 
outbursts. 

It is possible that this problem could be solved by adding an 
underlying nearly constant spectrum with a maximum around 
a frequency of 90 GHz, in association with a steeper optically 
thin spectral index for the individual outbursts. This would im- 
ply a higher value for the index s of the electron energy distri- 
bution, in better agreement with the value of 2.4 infeiTed from 
the measure of the millimetre-to-infrared slope of the flaring 
spectrum during the initial stage of the 1983 outburst (MG85). 
Such an underlying contribution from the inner jet to the mil- 
limetre and submillimetre emission of 3C 273 could possibly 
be related to the so called TZ component identified in the blue- 
bump of 3C 273 by Paltani et al. ( |1998| ). Indeed, this slowly 
varying component, which dominates the optical spectrum, can 
be interpreted as synchrotron emission and might reveal the 
blazar-like characteristics of 3C 273. 

As mentioned in Sect. ^ the obtained value of 1.6 for the 
index b suggests that the parallel and the perpendicular 
B± components of the magnetic field B are roughly equal. 
The simplest interpretation of this result is that the magnetic 
field B is turbulent behind the shock front, leading to a re- 
isotropization of the field. Another possibility, but which re- 
quires much more fine tuning, is that iJy is dominant in the 
underlying jet and that the shock amplifies just enough the B± 
component (see MG85) to end with B± « B\\. The obtained 
value of r ~ 0.8 suggests that the jet opens less with distance 
L along the jet axis than if it was conical. Such a behaviour 
was recently observed in M 87 (Virgo A, 3C 274) on an im- 
age of its very inner jet (Junor et al. 1999). A non-conical jet 
can arise due to an accelerating jet flow (Marscher 198C), but 
this interpretation would be in contradiction with the constant 
Doppler factor V assumed in this model. Another possibility 
is that there is an external pressure which confines the jet, so 
that the jet does not expand freely in two dimensions. Such a 
collimation could be of magnetic origin, either in the case of 



a magnetically self-confined jet (cf. Begelman et al. 1984) or 



in the case of the two-flow model (Sol et al. 1989), which as 



sumes that a highly relativistic electron-positron beam is con- 
fined within the magneto-hydrodynamic structure of a mildly 
relativistic electron-protonjet. 

The discussion above concerns the results obtained with a 
shock model in which we assume a constant Doppler factor T> 
by imposing a zero value for the index d. We obtain a com- 
parable of the fit using an alternative model in which we 
assume that the jet is conical and thus that the index r equals 
1 . For this model, we obtained very similar values for most pa- 
rameters given in Table ^ including the indices s and b, which 
therefore seem to be strongly constrained by the light-curve de- 
composition. What changes with this conical jet model are the 
values of the indices k and d. The value of k decreases to 2.72, 
which is now almost the value of fcad = 2(s+2) /3 ~ 2.70 cor- 
responding to an adiabatic jet flow and the best-fit value for d is 
0.10, which suggests either that the emitting region is slightly 
decelerating while it travels down the jet or that the jet is bend- 
ing away from the line-of-sight. Even if the shock front is not 



decelerating, a slight deceleration of the centre of the emitting 
region is a natural consequence of the continuous increase of its 
thickness x during the whole shock evolution. Therefore, even 
for a straight jet, this positive value of d might not be in con- 
tradiction with a value of k typical for an adiabatic jet flow. On 
the other hand, it is well established that the path projected on 
the sky followed by the VLBl components in the parsec-scale 
jet of 3C 273 is curved (e.g. Abraham et al. 1996). It is there- 



fore also possible that the slight decrease of the Doppler factor 

V suggested by the conical jet model is related to this overall 
bending of the jet. 

The fact that both the model with a constant Doppler factor 

V and the conical jet model give an equally good fit suggests 
that the light-curve decomposition cannot uniquely determine 
the values of the three parameters r, k and d. Indeed, by leaving 
all these parameters free to vary, we could not achieve a better 
fit than obtained with one of them fixed. It means that besides 
the two models discussed above, jet models with other combi- 
nations of the values of r, k and d cannot be excluded by our 
light-curve decomposition. 

6.2. On the peculiarities of individual outbursts 

Concerning the specificity of individual outbursts, we noted in 
Paper 1 that short- and long-lived outbursts in 3C 273 are usu- 
ally not peaking at the same frequency. Long-lived outbursts 
were found to peak at lower frequencies and we proposed that 
this relationship might be related to the distance down the jet 
at which the shock forms. The proper interpretation of the pe- 
culiarities of individual outbursts is now complicated by the 
fact that the outbursts are not self-similar anymore, because of 
the different behaviour of the two stage transitions 1 1 2 and 2 1 3 
when the physical quantities K, B and T) change (cf. Sect. 3.2). 
For instance, relatively high values of K^^ and low values of 
Son result in a shorter synchrotron stage. This effect is par- 
ticularly pronounced in the three successive outbursts between 
1982 and 1984 (cf. Table |). The apparent absence of a flat 
peaking stage in 3C 345 (Stevens et al. 1996) could therefore 



be due to a relatively low magnetic field B as compared to the 
factor K of the electron energy distribution. 

We note in Sect. |] that there is apparently no correlation 
from one outburst to the other between the two quantities A'on 
and i?on at the onset of the shock. This result suggests that 
differences from one outburst to the other are not primarily 
due to changes of the compression ratio of the shock, be- 
cause in this case we would expect correlated Kq^ and Son 
variations due to their similar dependence on 77 (cf. MG85). 
Changes of i^on which are independent of i?on could be pos- 
sible if the acceleration process of the electrons crossing the 
shock front was not adiabatic. In this case, the energy gain <^ 
of the electrons would not simply be related to the compres- 
sion ratio as ^ = t^'^/'^ (MG85), but might have different be- 
haviours from one outburst to the other. Since the normaliza- 
tion of the electron energy distribution is proportional to 
^^^^ (MG85), changes of the electron energy gain ^ for dif- 
ferent shocks in 3C 273 could lead to the observed changes 
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of KoTx- This interpretation is supported by the significant anti- 
correlation found between Ko^ and the Doppler factor I?on, be- 
cause a non-adiabatic shock wave converts bulk kinetic energy 



merical simulations of shocks in relativistic jets (e.g. Gomez et 
al. 1997) should lead in a near future to important progresses 



into internal energy (e.g. Begelman et al. 1984). Outbursts with 



a small value of V and a great value of K, would therefore be 
associated with more efficient shocks, for which a greater frac- 
tion of their bulk kinetic energy, as measured by T), would be 
converted into internal energy by increasing the average energy 
gain ^ of the electrons. 

7. Summary and conclusion 

This work presents a phenomenological model describing the 
evolution of synchrotron outbursts expected to be emitted by 
shock waves in a relativistic jet. This model, which is a gen- 
eralization of the original shock model of MG85, is then con- 
fronted to the very well sampled long-term submillimetre-to- 
radio Ught-curves of 3C 273. Many iterative fits are needed to 
adjust the more than 5000 observations with a series of sev- 
enteen successive outbursts defined by 85 parameters of the 
model. This fitting procedure allows us not only to define the 
average outburst's evolution in 3C 273, but also to study the pe- 
culiarities of each individual outburst during the last 20 years. 
The main results are the following: 

- The quiescent low frequency emission can be understood 
as being entirely due to the superimposition of slowly de- 
caying outbursts, which flared a few years before, and a 
constant contribution from the hot spot (3C 273A) of the 
outer jet. 

- The values of the indices describing how the physical quan- 
tities characterizing the jet evolve with its opening are 
found to be in general agreement with the simple jet model 
considered by MG85. 

- The best-fit values of these indices suggest however either 
that the jet opens slightly less than a conical jet if we as- 
sume that the emitting region has a constant bulk Doppler 
factor or, alternatively, that this Doppler factor decreases if 
we impose the jet to be conical. 

- In both cases, the magnetic field in the emitting region be- 
hind the shock front seems to be rather turbulent and the 
jet flow is found to be nearly adiabatic, especially for the 
conical jet model. 

- The peculiarities of individual outbursts can be understood 
as being due to changes at the onset of the shock of the 
magnetic field strength B, the normalization K of the elec- 
tron energy distribution and the Doppler factor T). 

- Shocks with a high value of K have usually a lower value 
of V. This anti-correlation might be related to the shock ef- 
ficiency to convert bulk kinetic energy into internal energy 
of the plasma behind the shock front. 

In the past, shock models were found to be difficult to test 
and constrain with total flux measurements. We show here that 
very strong observational constraints can be derived from long- 
term multi-wavelength monitoring campaigns. This new abil- 
ity, together with interferometric imaging techniques and nu- 



in the understanding of the physics involved in relativistic jets. 
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